3.16.36 \(\int \frac {1}{(d+e x)^3 (a^2+2 a b x+b^2 x^2)^3} \, dx\) [1536]

3.16.36.1 Optimal result
3.16.36.2 Mathematica [A] (verified)
3.16.36.3 Rubi [A] (verified)
3.16.36.4 Maple [A] (verified)
3.16.36.5 Fricas [B] (verification not implemented)
3.16.36.6 Sympy [B] (verification not implemented)
3.16.36.7 Maxima [B] (verification not implemented)
3.16.36.8 Giac [B] (verification not implemented)
3.16.36.9 Mupad [B] (verification not implemented)

3.16.36.1 Optimal result

Integrand size = 26, antiderivative size = 220 \[ \int \frac {1}{(d+e x)^3 \left (a^2+2 a b x+b^2 x^2\right )^3} \, dx=-\frac {b^2}{5 (b d-a e)^3 (a+b x)^5}+\frac {3 b^2 e}{4 (b d-a e)^4 (a+b x)^4}-\frac {2 b^2 e^2}{(b d-a e)^5 (a+b x)^3}+\frac {5 b^2 e^3}{(b d-a e)^6 (a+b x)^2}-\frac {15 b^2 e^4}{(b d-a e)^7 (a+b x)}-\frac {e^5}{2 (b d-a e)^6 (d+e x)^2}-\frac {6 b e^5}{(b d-a e)^7 (d+e x)}-\frac {21 b^2 e^5 \log (a+b x)}{(b d-a e)^8}+\frac {21 b^2 e^5 \log (d+e x)}{(b d-a e)^8} \]

output
-1/5*b^2/(-a*e+b*d)^3/(b*x+a)^5+3/4*b^2*e/(-a*e+b*d)^4/(b*x+a)^4-2*b^2*e^2 
/(-a*e+b*d)^5/(b*x+a)^3+5*b^2*e^3/(-a*e+b*d)^6/(b*x+a)^2-15*b^2*e^4/(-a*e+ 
b*d)^7/(b*x+a)-1/2*e^5/(-a*e+b*d)^6/(e*x+d)^2-6*b*e^5/(-a*e+b*d)^7/(e*x+d) 
-21*b^2*e^5*ln(b*x+a)/(-a*e+b*d)^8+21*b^2*e^5*ln(e*x+d)/(-a*e+b*d)^8
 
3.16.36.2 Mathematica [A] (verified)

Time = 0.10 (sec) , antiderivative size = 204, normalized size of antiderivative = 0.93 \[ \int \frac {1}{(d+e x)^3 \left (a^2+2 a b x+b^2 x^2\right )^3} \, dx=-\frac {\frac {4 b^2 (b d-a e)^5}{(a+b x)^5}-\frac {15 b^2 e (b d-a e)^4}{(a+b x)^4}+\frac {40 b^2 e^2 (b d-a e)^3}{(a+b x)^3}-\frac {100 b^2 e^3 (b d-a e)^2}{(a+b x)^2}+\frac {300 b^2 e^4 (b d-a e)}{a+b x}+\frac {10 e^5 (b d-a e)^2}{(d+e x)^2}+\frac {120 b e^5 (b d-a e)}{d+e x}+420 b^2 e^5 \log (a+b x)-420 b^2 e^5 \log (d+e x)}{20 (b d-a e)^8} \]

input
Integrate[1/((d + e*x)^3*(a^2 + 2*a*b*x + b^2*x^2)^3),x]
 
output
-1/20*((4*b^2*(b*d - a*e)^5)/(a + b*x)^5 - (15*b^2*e*(b*d - a*e)^4)/(a + b 
*x)^4 + (40*b^2*e^2*(b*d - a*e)^3)/(a + b*x)^3 - (100*b^2*e^3*(b*d - a*e)^ 
2)/(a + b*x)^2 + (300*b^2*e^4*(b*d - a*e))/(a + b*x) + (10*e^5*(b*d - a*e) 
^2)/(d + e*x)^2 + (120*b*e^5*(b*d - a*e))/(d + e*x) + 420*b^2*e^5*Log[a + 
b*x] - 420*b^2*e^5*Log[d + e*x])/(b*d - a*e)^8
 
3.16.36.3 Rubi [A] (verified)

Time = 0.48 (sec) , antiderivative size = 220, normalized size of antiderivative = 1.00, number of steps used = 4, number of rules used = 4, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.154, Rules used = {1098, 27, 54, 2009}

Below are the steps used by Rubi to obtain the solution. The rule number used for the transformation is given above next to the arrow. The rules definitions used are listed below.

\(\displaystyle \int \frac {1}{\left (a^2+2 a b x+b^2 x^2\right )^3 (d+e x)^3} \, dx\)

\(\Big \downarrow \) 1098

\(\displaystyle b^6 \int \frac {1}{b^6 (a+b x)^6 (d+e x)^3}dx\)

\(\Big \downarrow \) 27

\(\displaystyle \int \frac {1}{(a+b x)^6 (d+e x)^3}dx\)

\(\Big \downarrow \) 54

\(\displaystyle \int \left (-\frac {21 b^3 e^5}{(a+b x) (b d-a e)^8}+\frac {15 b^3 e^4}{(a+b x)^2 (b d-a e)^7}-\frac {10 b^3 e^3}{(a+b x)^3 (b d-a e)^6}+\frac {6 b^3 e^2}{(a+b x)^4 (b d-a e)^5}-\frac {3 b^3 e}{(a+b x)^5 (b d-a e)^4}+\frac {b^3}{(a+b x)^6 (b d-a e)^3}+\frac {21 b^2 e^6}{(d+e x) (b d-a e)^8}+\frac {6 b e^6}{(d+e x)^2 (b d-a e)^7}+\frac {e^6}{(d+e x)^3 (b d-a e)^6}\right )dx\)

\(\Big \downarrow \) 2009

\(\displaystyle -\frac {21 b^2 e^5 \log (a+b x)}{(b d-a e)^8}+\frac {21 b^2 e^5 \log (d+e x)}{(b d-a e)^8}-\frac {15 b^2 e^4}{(a+b x) (b d-a e)^7}+\frac {5 b^2 e^3}{(a+b x)^2 (b d-a e)^6}-\frac {2 b^2 e^2}{(a+b x)^3 (b d-a e)^5}+\frac {3 b^2 e}{4 (a+b x)^4 (b d-a e)^4}-\frac {b^2}{5 (a+b x)^5 (b d-a e)^3}-\frac {6 b e^5}{(d+e x) (b d-a e)^7}-\frac {e^5}{2 (d+e x)^2 (b d-a e)^6}\)

input
Int[1/((d + e*x)^3*(a^2 + 2*a*b*x + b^2*x^2)^3),x]
 
output
-1/5*b^2/((b*d - a*e)^3*(a + b*x)^5) + (3*b^2*e)/(4*(b*d - a*e)^4*(a + b*x 
)^4) - (2*b^2*e^2)/((b*d - a*e)^5*(a + b*x)^3) + (5*b^2*e^3)/((b*d - a*e)^ 
6*(a + b*x)^2) - (15*b^2*e^4)/((b*d - a*e)^7*(a + b*x)) - e^5/(2*(b*d - a* 
e)^6*(d + e*x)^2) - (6*b*e^5)/((b*d - a*e)^7*(d + e*x)) - (21*b^2*e^5*Log[ 
a + b*x])/(b*d - a*e)^8 + (21*b^2*e^5*Log[d + e*x])/(b*d - a*e)^8
 

3.16.36.3.1 Defintions of rubi rules used

rule 27
Int[(a_)*(Fx_), x_Symbol] :> Simp[a   Int[Fx, x], x] /; FreeQ[a, x] &&  !Ma 
tchQ[Fx, (b_)*(Gx_) /; FreeQ[b, x]]
 

rule 54
Int[((a_) + (b_.)*(x_))^(m_)*((c_.) + (d_.)*(x_))^(n_.), x_Symbol] :> Int[E 
xpandIntegrand[(a + b*x)^m*(c + d*x)^n, x], x] /; FreeQ[{a, b, c, d}, x] && 
 ILtQ[m, 0] && IntegerQ[n] &&  !(IGtQ[n, 0] && LtQ[m + n + 2, 0])
 

rule 1098
Int[((d_.) + (e_.)*(x_))^(m_.)*((a_) + (b_.)*(x_) + (c_.)*(x_)^2)^(p_.), x_ 
Symbol] :> Simp[1/c^p   Int[(d + e*x)^m*(b/2 + c*x)^(2*p), x], x] /; FreeQ[ 
{a, b, c, d, e, m}, x] && EqQ[b^2 - 4*a*c, 0] && IntegerQ[p]
 

rule 2009
Int[u_, x_Symbol] :> Simp[IntSum[u, x], x] /; SumQ[u]
 
3.16.36.4 Maple [A] (verified)

Time = 2.36 (sec) , antiderivative size = 215, normalized size of antiderivative = 0.98

method result size
default \(\frac {b^{2}}{5 \left (a e -b d \right )^{3} \left (b x +a \right )^{5}}-\frac {21 b^{2} e^{5} \ln \left (b x +a \right )}{\left (a e -b d \right )^{8}}+\frac {15 b^{2} e^{4}}{\left (a e -b d \right )^{7} \left (b x +a \right )}+\frac {5 b^{2} e^{3}}{\left (a e -b d \right )^{6} \left (b x +a \right )^{2}}+\frac {2 b^{2} e^{2}}{\left (a e -b d \right )^{5} \left (b x +a \right )^{3}}+\frac {3 b^{2} e}{4 \left (a e -b d \right )^{4} \left (b x +a \right )^{4}}-\frac {e^{5}}{2 \left (a e -b d \right )^{6} \left (e x +d \right )^{2}}+\frac {21 b^{2} e^{5} \ln \left (e x +d \right )}{\left (a e -b d \right )^{8}}+\frac {6 e^{5} b}{\left (a e -b d \right )^{7} \left (e x +d \right )}\) \(215\)
risch \(\text {Expression too large to display}\) \(1264\)
norman \(\text {Expression too large to display}\) \(1314\)
parallelrisch \(\text {Expression too large to display}\) \(1360\)

input
int(1/(e*x+d)^3/(b^2*x^2+2*a*b*x+a^2)^3,x,method=_RETURNVERBOSE)
 
output
1/5*b^2/(a*e-b*d)^3/(b*x+a)^5-21*b^2/(a*e-b*d)^8*e^5*ln(b*x+a)+15*b^2/(a*e 
-b*d)^7*e^4/(b*x+a)+5*b^2/(a*e-b*d)^6*e^3/(b*x+a)^2+2*b^2/(a*e-b*d)^5*e^2/ 
(b*x+a)^3+3/4*b^2/(a*e-b*d)^4*e/(b*x+a)^4-1/2*e^5/(a*e-b*d)^6/(e*x+d)^2+21 
*b^2/(a*e-b*d)^8*e^5*ln(e*x+d)+6*e^5/(a*e-b*d)^7*b/(e*x+d)
 
3.16.36.5 Fricas [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 2005 vs. \(2 (214) = 428\).

Time = 0.38 (sec) , antiderivative size = 2005, normalized size of antiderivative = 9.11 \[ \int \frac {1}{(d+e x)^3 \left (a^2+2 a b x+b^2 x^2\right )^3} \, dx=\text {Too large to display} \]

input
integrate(1/(e*x+d)^3/(b^2*x^2+2*a*b*x+a^2)^3,x, algorithm="fricas")
 
output
-1/20*(4*b^7*d^7 - 35*a*b^6*d^6*e + 140*a^2*b^5*d^5*e^2 - 350*a^3*b^4*d^4* 
e^3 + 700*a^4*b^3*d^3*e^4 - 329*a^5*b^2*d^2*e^5 - 140*a^6*b*d*e^6 + 10*a^7 
*e^7 + 420*(b^7*d*e^6 - a*b^6*e^7)*x^6 + 630*(b^7*d^2*e^5 + 2*a*b^6*d*e^6 
- 3*a^2*b^5*e^7)*x^5 + 70*(2*b^7*d^3*e^4 + 39*a*b^6*d^2*e^5 + 6*a^2*b^5*d* 
e^6 - 47*a^3*b^4*e^7)*x^4 - 35*(b^7*d^4*e^3 - 20*a*b^6*d^3*e^4 - 126*a^2*b 
^5*d^2*e^5 + 68*a^3*b^4*d*e^6 + 77*a^4*b^3*e^7)*x^3 + 7*(2*b^7*d^5*e^2 - 2 
5*a*b^6*d^4*e^3 + 200*a^2*b^5*d^3*e^4 + 430*a^3*b^4*d^2*e^5 - 470*a^4*b^3* 
d*e^6 - 137*a^5*b^2*e^7)*x^2 - 7*(b^7*d^6*e - 10*a*b^6*d^5*e^2 + 50*a^2*b^ 
5*d^4*e^3 - 200*a^3*b^4*d^3*e^4 - 65*a^4*b^3*d^2*e^5 + 214*a^5*b^2*d*e^6 + 
 10*a^6*b*e^7)*x + 420*(b^7*e^7*x^7 + a^5*b^2*d^2*e^5 + (2*b^7*d*e^6 + 5*a 
*b^6*e^7)*x^6 + (b^7*d^2*e^5 + 10*a*b^6*d*e^6 + 10*a^2*b^5*e^7)*x^5 + 5*(a 
*b^6*d^2*e^5 + 4*a^2*b^5*d*e^6 + 2*a^3*b^4*e^7)*x^4 + 5*(2*a^2*b^5*d^2*e^5 
 + 4*a^3*b^4*d*e^6 + a^4*b^3*e^7)*x^3 + (10*a^3*b^4*d^2*e^5 + 10*a^4*b^3*d 
*e^6 + a^5*b^2*e^7)*x^2 + (5*a^4*b^3*d^2*e^5 + 2*a^5*b^2*d*e^6)*x)*log(b*x 
 + a) - 420*(b^7*e^7*x^7 + a^5*b^2*d^2*e^5 + (2*b^7*d*e^6 + 5*a*b^6*e^7)*x 
^6 + (b^7*d^2*e^5 + 10*a*b^6*d*e^6 + 10*a^2*b^5*e^7)*x^5 + 5*(a*b^6*d^2*e^ 
5 + 4*a^2*b^5*d*e^6 + 2*a^3*b^4*e^7)*x^4 + 5*(2*a^2*b^5*d^2*e^5 + 4*a^3*b^ 
4*d*e^6 + a^4*b^3*e^7)*x^3 + (10*a^3*b^4*d^2*e^5 + 10*a^4*b^3*d*e^6 + a^5* 
b^2*e^7)*x^2 + (5*a^4*b^3*d^2*e^5 + 2*a^5*b^2*d*e^6)*x)*log(e*x + d))/(a^5 
*b^8*d^10 - 8*a^6*b^7*d^9*e + 28*a^7*b^6*d^8*e^2 - 56*a^8*b^5*d^7*e^3 +...
 
3.16.36.6 Sympy [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 1974 vs. \(2 (201) = 402\).

Time = 39.81 (sec) , antiderivative size = 1974, normalized size of antiderivative = 8.97 \[ \int \frac {1}{(d+e x)^3 \left (a^2+2 a b x+b^2 x^2\right )^3} \, dx=\text {Too large to display} \]

input
integrate(1/(e*x+d)**3/(b**2*x**2+2*a*b*x+a**2)**3,x)
 
output
21*b**2*e**5*log(x + (-21*a**9*b**2*e**14/(a*e - b*d)**8 + 189*a**8*b**3*d 
*e**13/(a*e - b*d)**8 - 756*a**7*b**4*d**2*e**12/(a*e - b*d)**8 + 1764*a** 
6*b**5*d**3*e**11/(a*e - b*d)**8 - 2646*a**5*b**6*d**4*e**10/(a*e - b*d)** 
8 + 2646*a**4*b**7*d**5*e**9/(a*e - b*d)**8 - 1764*a**3*b**8*d**6*e**8/(a* 
e - b*d)**8 + 756*a**2*b**9*d**7*e**7/(a*e - b*d)**8 - 189*a*b**10*d**8*e* 
*6/(a*e - b*d)**8 + 21*a*b**2*e**6 + 21*b**11*d**9*e**5/(a*e - b*d)**8 + 2 
1*b**3*d*e**5)/(42*b**3*e**6))/(a*e - b*d)**8 - 21*b**2*e**5*log(x + (21*a 
**9*b**2*e**14/(a*e - b*d)**8 - 189*a**8*b**3*d*e**13/(a*e - b*d)**8 + 756 
*a**7*b**4*d**2*e**12/(a*e - b*d)**8 - 1764*a**6*b**5*d**3*e**11/(a*e - b* 
d)**8 + 2646*a**5*b**6*d**4*e**10/(a*e - b*d)**8 - 2646*a**4*b**7*d**5*e** 
9/(a*e - b*d)**8 + 1764*a**3*b**8*d**6*e**8/(a*e - b*d)**8 - 756*a**2*b**9 
*d**7*e**7/(a*e - b*d)**8 + 189*a*b**10*d**8*e**6/(a*e - b*d)**8 + 21*a*b* 
*2*e**6 - 21*b**11*d**9*e**5/(a*e - b*d)**8 + 21*b**3*d*e**5)/(42*b**3*e** 
6))/(a*e - b*d)**8 + (-10*a**6*e**6 + 130*a**5*b*d*e**5 + 459*a**4*b**2*d* 
*2*e**4 - 241*a**3*b**3*d**3*e**3 + 109*a**2*b**4*d**4*e**2 - 31*a*b**5*d* 
*5*e + 4*b**6*d**6 + 420*b**6*e**6*x**6 + x**5*(1890*a*b**5*e**6 + 630*b** 
6*d*e**5) + x**4*(3290*a**2*b**4*e**6 + 2870*a*b**5*d*e**5 + 140*b**6*d**2 
*e**4) + x**3*(2695*a**3*b**3*e**6 + 5075*a**2*b**4*d*e**5 + 665*a*b**5*d* 
*2*e**4 - 35*b**6*d**3*e**3) + x**2*(959*a**4*b**2*e**6 + 4249*a**3*b**3*d 
*e**5 + 1239*a**2*b**4*d**2*e**4 - 161*a*b**5*d**3*e**3 + 14*b**6*d**4*...
 
3.16.36.7 Maxima [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 1558 vs. \(2 (214) = 428\).

Time = 0.30 (sec) , antiderivative size = 1558, normalized size of antiderivative = 7.08 \[ \int \frac {1}{(d+e x)^3 \left (a^2+2 a b x+b^2 x^2\right )^3} \, dx=\text {Too large to display} \]

input
integrate(1/(e*x+d)^3/(b^2*x^2+2*a*b*x+a^2)^3,x, algorithm="maxima")
 
output
-21*b^2*e^5*log(b*x + a)/(b^8*d^8 - 8*a*b^7*d^7*e + 28*a^2*b^6*d^6*e^2 - 5 
6*a^3*b^5*d^5*e^3 + 70*a^4*b^4*d^4*e^4 - 56*a^5*b^3*d^3*e^5 + 28*a^6*b^2*d 
^2*e^6 - 8*a^7*b*d*e^7 + a^8*e^8) + 21*b^2*e^5*log(e*x + d)/(b^8*d^8 - 8*a 
*b^7*d^7*e + 28*a^2*b^6*d^6*e^2 - 56*a^3*b^5*d^5*e^3 + 70*a^4*b^4*d^4*e^4 
- 56*a^5*b^3*d^3*e^5 + 28*a^6*b^2*d^2*e^6 - 8*a^7*b*d*e^7 + a^8*e^8) - 1/2 
0*(420*b^6*e^6*x^6 + 4*b^6*d^6 - 31*a*b^5*d^5*e + 109*a^2*b^4*d^4*e^2 - 24 
1*a^3*b^3*d^3*e^3 + 459*a^4*b^2*d^2*e^4 + 130*a^5*b*d*e^5 - 10*a^6*e^6 + 6 
30*(b^6*d*e^5 + 3*a*b^5*e^6)*x^5 + 70*(2*b^6*d^2*e^4 + 41*a*b^5*d*e^5 + 47 
*a^2*b^4*e^6)*x^4 - 35*(b^6*d^3*e^3 - 19*a*b^5*d^2*e^4 - 145*a^2*b^4*d*e^5 
 - 77*a^3*b^3*e^6)*x^3 + 7*(2*b^6*d^4*e^2 - 23*a*b^5*d^3*e^3 + 177*a^2*b^4 
*d^2*e^4 + 607*a^3*b^3*d*e^5 + 137*a^4*b^2*e^6)*x^2 - 7*(b^6*d^5*e - 9*a*b 
^5*d^4*e^2 + 41*a^2*b^4*d^3*e^3 - 159*a^3*b^3*d^2*e^4 - 224*a^4*b^2*d*e^5 
- 10*a^5*b*e^6)*x)/(a^5*b^7*d^9 - 7*a^6*b^6*d^8*e + 21*a^7*b^5*d^7*e^2 - 3 
5*a^8*b^4*d^6*e^3 + 35*a^9*b^3*d^5*e^4 - 21*a^10*b^2*d^4*e^5 + 7*a^11*b*d^ 
3*e^6 - a^12*d^2*e^7 + (b^12*d^7*e^2 - 7*a*b^11*d^6*e^3 + 21*a^2*b^10*d^5* 
e^4 - 35*a^3*b^9*d^4*e^5 + 35*a^4*b^8*d^3*e^6 - 21*a^5*b^7*d^2*e^7 + 7*a^6 
*b^6*d*e^8 - a^7*b^5*e^9)*x^7 + (2*b^12*d^8*e - 9*a*b^11*d^7*e^2 + 7*a^2*b 
^10*d^6*e^3 + 35*a^3*b^9*d^5*e^4 - 105*a^4*b^8*d^4*e^5 + 133*a^5*b^7*d^3*e 
^6 - 91*a^6*b^6*d^2*e^7 + 33*a^7*b^5*d*e^8 - 5*a^8*b^4*e^9)*x^6 + (b^12*d^ 
9 + 3*a*b^11*d^8*e - 39*a^2*b^10*d^7*e^2 + 105*a^3*b^9*d^6*e^3 - 105*a^...
 
3.16.36.8 Giac [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 714 vs. \(2 (214) = 428\).

Time = 0.27 (sec) , antiderivative size = 714, normalized size of antiderivative = 3.25 \[ \int \frac {1}{(d+e x)^3 \left (a^2+2 a b x+b^2 x^2\right )^3} \, dx=-\frac {21 \, b^{3} e^{5} \log \left ({\left | b x + a \right |}\right )}{b^{9} d^{8} - 8 \, a b^{8} d^{7} e + 28 \, a^{2} b^{7} d^{6} e^{2} - 56 \, a^{3} b^{6} d^{5} e^{3} + 70 \, a^{4} b^{5} d^{4} e^{4} - 56 \, a^{5} b^{4} d^{3} e^{5} + 28 \, a^{6} b^{3} d^{2} e^{6} - 8 \, a^{7} b^{2} d e^{7} + a^{8} b e^{8}} + \frac {21 \, b^{2} e^{6} \log \left ({\left | e x + d \right |}\right )}{b^{8} d^{8} e - 8 \, a b^{7} d^{7} e^{2} + 28 \, a^{2} b^{6} d^{6} e^{3} - 56 \, a^{3} b^{5} d^{5} e^{4} + 70 \, a^{4} b^{4} d^{4} e^{5} - 56 \, a^{5} b^{3} d^{3} e^{6} + 28 \, a^{6} b^{2} d^{2} e^{7} - 8 \, a^{7} b d e^{8} + a^{8} e^{9}} - \frac {4 \, b^{7} d^{7} - 35 \, a b^{6} d^{6} e + 140 \, a^{2} b^{5} d^{5} e^{2} - 350 \, a^{3} b^{4} d^{4} e^{3} + 700 \, a^{4} b^{3} d^{3} e^{4} - 329 \, a^{5} b^{2} d^{2} e^{5} - 140 \, a^{6} b d e^{6} + 10 \, a^{7} e^{7} + 420 \, {\left (b^{7} d e^{6} - a b^{6} e^{7}\right )} x^{6} + 630 \, {\left (b^{7} d^{2} e^{5} + 2 \, a b^{6} d e^{6} - 3 \, a^{2} b^{5} e^{7}\right )} x^{5} + 70 \, {\left (2 \, b^{7} d^{3} e^{4} + 39 \, a b^{6} d^{2} e^{5} + 6 \, a^{2} b^{5} d e^{6} - 47 \, a^{3} b^{4} e^{7}\right )} x^{4} - 35 \, {\left (b^{7} d^{4} e^{3} - 20 \, a b^{6} d^{3} e^{4} - 126 \, a^{2} b^{5} d^{2} e^{5} + 68 \, a^{3} b^{4} d e^{6} + 77 \, a^{4} b^{3} e^{7}\right )} x^{3} + 7 \, {\left (2 \, b^{7} d^{5} e^{2} - 25 \, a b^{6} d^{4} e^{3} + 200 \, a^{2} b^{5} d^{3} e^{4} + 430 \, a^{3} b^{4} d^{2} e^{5} - 470 \, a^{4} b^{3} d e^{6} - 137 \, a^{5} b^{2} e^{7}\right )} x^{2} - 7 \, {\left (b^{7} d^{6} e - 10 \, a b^{6} d^{5} e^{2} + 50 \, a^{2} b^{5} d^{4} e^{3} - 200 \, a^{3} b^{4} d^{3} e^{4} - 65 \, a^{4} b^{3} d^{2} e^{5} + 214 \, a^{5} b^{2} d e^{6} + 10 \, a^{6} b e^{7}\right )} x}{20 \, {\left (b d - a e\right )}^{8} {\left (b x + a\right )}^{5} {\left (e x + d\right )}^{2}} \]

input
integrate(1/(e*x+d)^3/(b^2*x^2+2*a*b*x+a^2)^3,x, algorithm="giac")
 
output
-21*b^3*e^5*log(abs(b*x + a))/(b^9*d^8 - 8*a*b^8*d^7*e + 28*a^2*b^7*d^6*e^ 
2 - 56*a^3*b^6*d^5*e^3 + 70*a^4*b^5*d^4*e^4 - 56*a^5*b^4*d^3*e^5 + 28*a^6* 
b^3*d^2*e^6 - 8*a^7*b^2*d*e^7 + a^8*b*e^8) + 21*b^2*e^6*log(abs(e*x + d))/ 
(b^8*d^8*e - 8*a*b^7*d^7*e^2 + 28*a^2*b^6*d^6*e^3 - 56*a^3*b^5*d^5*e^4 + 7 
0*a^4*b^4*d^4*e^5 - 56*a^5*b^3*d^3*e^6 + 28*a^6*b^2*d^2*e^7 - 8*a^7*b*d*e^ 
8 + a^8*e^9) - 1/20*(4*b^7*d^7 - 35*a*b^6*d^6*e + 140*a^2*b^5*d^5*e^2 - 35 
0*a^3*b^4*d^4*e^3 + 700*a^4*b^3*d^3*e^4 - 329*a^5*b^2*d^2*e^5 - 140*a^6*b* 
d*e^6 + 10*a^7*e^7 + 420*(b^7*d*e^6 - a*b^6*e^7)*x^6 + 630*(b^7*d^2*e^5 + 
2*a*b^6*d*e^6 - 3*a^2*b^5*e^7)*x^5 + 70*(2*b^7*d^3*e^4 + 39*a*b^6*d^2*e^5 
+ 6*a^2*b^5*d*e^6 - 47*a^3*b^4*e^7)*x^4 - 35*(b^7*d^4*e^3 - 20*a*b^6*d^3*e 
^4 - 126*a^2*b^5*d^2*e^5 + 68*a^3*b^4*d*e^6 + 77*a^4*b^3*e^7)*x^3 + 7*(2*b 
^7*d^5*e^2 - 25*a*b^6*d^4*e^3 + 200*a^2*b^5*d^3*e^4 + 430*a^3*b^4*d^2*e^5 
- 470*a^4*b^3*d*e^6 - 137*a^5*b^2*e^7)*x^2 - 7*(b^7*d^6*e - 10*a*b^6*d^5*e 
^2 + 50*a^2*b^5*d^4*e^3 - 200*a^3*b^4*d^3*e^4 - 65*a^4*b^3*d^2*e^5 + 214*a 
^5*b^2*d*e^6 + 10*a^6*b*e^7)*x)/((b*d - a*e)^8*(b*x + a)^5*(e*x + d)^2)
 
3.16.36.9 Mupad [B] (verification not implemented)

Time = 10.58 (sec) , antiderivative size = 1427, normalized size of antiderivative = 6.49 \[ \int \frac {1}{(d+e x)^3 \left (a^2+2 a b x+b^2 x^2\right )^3} \, dx=\frac {\frac {-10\,a^6\,e^6+130\,a^5\,b\,d\,e^5+459\,a^4\,b^2\,d^2\,e^4-241\,a^3\,b^3\,d^3\,e^3+109\,a^2\,b^4\,d^4\,e^2-31\,a\,b^5\,d^5\,e+4\,b^6\,d^6}{20\,\left (a^7\,e^7-7\,a^6\,b\,d\,e^6+21\,a^5\,b^2\,d^2\,e^5-35\,a^4\,b^3\,d^3\,e^4+35\,a^3\,b^4\,d^4\,e^3-21\,a^2\,b^5\,d^5\,e^2+7\,a\,b^6\,d^6\,e-b^7\,d^7\right )}+\frac {7\,e^3\,x^3\,\left (77\,a^3\,b^3\,e^3+145\,a^2\,b^4\,d\,e^2+19\,a\,b^5\,d^2\,e-b^6\,d^3\right )}{4\,\left (a^7\,e^7-7\,a^6\,b\,d\,e^6+21\,a^5\,b^2\,d^2\,e^5-35\,a^4\,b^3\,d^3\,e^4+35\,a^3\,b^4\,d^4\,e^3-21\,a^2\,b^5\,d^5\,e^2+7\,a\,b^6\,d^6\,e-b^7\,d^7\right )}+\frac {21\,b^6\,e^6\,x^6}{a^7\,e^7-7\,a^6\,b\,d\,e^6+21\,a^5\,b^2\,d^2\,e^5-35\,a^4\,b^3\,d^3\,e^4+35\,a^3\,b^4\,d^4\,e^3-21\,a^2\,b^5\,d^5\,e^2+7\,a\,b^6\,d^6\,e-b^7\,d^7}+\frac {7\,e^2\,x^2\,\left (137\,a^4\,b^2\,e^4+607\,a^3\,b^3\,d\,e^3+177\,a^2\,b^4\,d^2\,e^2-23\,a\,b^5\,d^3\,e+2\,b^6\,d^4\right )}{20\,\left (a^7\,e^7-7\,a^6\,b\,d\,e^6+21\,a^5\,b^2\,d^2\,e^5-35\,a^4\,b^3\,d^3\,e^4+35\,a^3\,b^4\,d^4\,e^3-21\,a^2\,b^5\,d^5\,e^2+7\,a\,b^6\,d^6\,e-b^7\,d^7\right )}+\frac {7\,e^4\,x^4\,\left (47\,a^2\,b^4\,e^2+41\,a\,b^5\,d\,e+2\,b^6\,d^2\right )}{2\,\left (a^7\,e^7-7\,a^6\,b\,d\,e^6+21\,a^5\,b^2\,d^2\,e^5-35\,a^4\,b^3\,d^3\,e^4+35\,a^3\,b^4\,d^4\,e^3-21\,a^2\,b^5\,d^5\,e^2+7\,a\,b^6\,d^6\,e-b^7\,d^7\right )}+\frac {7\,e\,x\,\left (10\,a^5\,b\,e^5+224\,a^4\,b^2\,d\,e^4+159\,a^3\,b^3\,d^2\,e^3-41\,a^2\,b^4\,d^3\,e^2+9\,a\,b^5\,d^4\,e-b^6\,d^5\right )}{20\,\left (a^7\,e^7-7\,a^6\,b\,d\,e^6+21\,a^5\,b^2\,d^2\,e^5-35\,a^4\,b^3\,d^3\,e^4+35\,a^3\,b^4\,d^4\,e^3-21\,a^2\,b^5\,d^5\,e^2+7\,a\,b^6\,d^6\,e-b^7\,d^7\right )}+\frac {63\,b\,e^4\,x^5\,\left (d\,b^5\,e+3\,a\,b^4\,e^2\right )}{2\,\left (a^7\,e^7-7\,a^6\,b\,d\,e^6+21\,a^5\,b^2\,d^2\,e^5-35\,a^4\,b^3\,d^3\,e^4+35\,a^3\,b^4\,d^4\,e^3-21\,a^2\,b^5\,d^5\,e^2+7\,a\,b^6\,d^6\,e-b^7\,d^7\right )}}{x^3\,\left (5\,a^4\,b\,e^2+20\,a^3\,b^2\,d\,e+10\,a^2\,b^3\,d^2\right )+x^4\,\left (10\,a^3\,b^2\,e^2+20\,a^2\,b^3\,d\,e+5\,a\,b^4\,d^2\right )+x\,\left (2\,e\,a^5\,d+5\,b\,a^4\,d^2\right )+x^2\,\left (a^5\,e^2+10\,a^4\,b\,d\,e+10\,a^3\,b^2\,d^2\right )+x^5\,\left (10\,a^2\,b^3\,e^2+10\,a\,b^4\,d\,e+b^5\,d^2\right )+x^6\,\left (2\,d\,b^5\,e+5\,a\,b^4\,e^2\right )+a^5\,d^2+b^5\,e^2\,x^7}-\frac {42\,b^2\,e^5\,\mathrm {atanh}\left (\frac {a^8\,e^8-6\,a^7\,b\,d\,e^7+14\,a^6\,b^2\,d^2\,e^6-14\,a^5\,b^3\,d^3\,e^5+14\,a^3\,b^5\,d^5\,e^3-14\,a^2\,b^6\,d^6\,e^2+6\,a\,b^7\,d^7\,e-b^8\,d^8}{{\left (a\,e-b\,d\right )}^8}+\frac {2\,b\,e\,x\,\left (a^7\,e^7-7\,a^6\,b\,d\,e^6+21\,a^5\,b^2\,d^2\,e^5-35\,a^4\,b^3\,d^3\,e^4+35\,a^3\,b^4\,d^4\,e^3-21\,a^2\,b^5\,d^5\,e^2+7\,a\,b^6\,d^6\,e-b^7\,d^7\right )}{{\left (a\,e-b\,d\right )}^8}\right )}{{\left (a\,e-b\,d\right )}^8} \]

input
int(1/((d + e*x)^3*(a^2 + b^2*x^2 + 2*a*b*x)^3),x)
 
output
((4*b^6*d^6 - 10*a^6*e^6 + 109*a^2*b^4*d^4*e^2 - 241*a^3*b^3*d^3*e^3 + 459 
*a^4*b^2*d^2*e^4 - 31*a*b^5*d^5*e + 130*a^5*b*d*e^5)/(20*(a^7*e^7 - b^7*d^ 
7 - 21*a^2*b^5*d^5*e^2 + 35*a^3*b^4*d^4*e^3 - 35*a^4*b^3*d^3*e^4 + 21*a^5* 
b^2*d^2*e^5 + 7*a*b^6*d^6*e - 7*a^6*b*d*e^6)) + (7*e^3*x^3*(77*a^3*b^3*e^3 
 - b^6*d^3 + 145*a^2*b^4*d*e^2 + 19*a*b^5*d^2*e))/(4*(a^7*e^7 - b^7*d^7 - 
21*a^2*b^5*d^5*e^2 + 35*a^3*b^4*d^4*e^3 - 35*a^4*b^3*d^3*e^4 + 21*a^5*b^2* 
d^2*e^5 + 7*a*b^6*d^6*e - 7*a^6*b*d*e^6)) + (21*b^6*e^6*x^6)/(a^7*e^7 - b^ 
7*d^7 - 21*a^2*b^5*d^5*e^2 + 35*a^3*b^4*d^4*e^3 - 35*a^4*b^3*d^3*e^4 + 21* 
a^5*b^2*d^2*e^5 + 7*a*b^6*d^6*e - 7*a^6*b*d*e^6) + (7*e^2*x^2*(2*b^6*d^4 + 
 137*a^4*b^2*e^4 + 607*a^3*b^3*d*e^3 + 177*a^2*b^4*d^2*e^2 - 23*a*b^5*d^3* 
e))/(20*(a^7*e^7 - b^7*d^7 - 21*a^2*b^5*d^5*e^2 + 35*a^3*b^4*d^4*e^3 - 35* 
a^4*b^3*d^3*e^4 + 21*a^5*b^2*d^2*e^5 + 7*a*b^6*d^6*e - 7*a^6*b*d*e^6)) + ( 
7*e^4*x^4*(2*b^6*d^2 + 47*a^2*b^4*e^2 + 41*a*b^5*d*e))/(2*(a^7*e^7 - b^7*d 
^7 - 21*a^2*b^5*d^5*e^2 + 35*a^3*b^4*d^4*e^3 - 35*a^4*b^3*d^3*e^4 + 21*a^5 
*b^2*d^2*e^5 + 7*a*b^6*d^6*e - 7*a^6*b*d*e^6)) + (7*e*x*(10*a^5*b*e^5 - b^ 
6*d^5 + 224*a^4*b^2*d*e^4 - 41*a^2*b^4*d^3*e^2 + 159*a^3*b^3*d^2*e^3 + 9*a 
*b^5*d^4*e))/(20*(a^7*e^7 - b^7*d^7 - 21*a^2*b^5*d^5*e^2 + 35*a^3*b^4*d^4* 
e^3 - 35*a^4*b^3*d^3*e^4 + 21*a^5*b^2*d^2*e^5 + 7*a*b^6*d^6*e - 7*a^6*b*d* 
e^6)) + (63*b*e^4*x^5*(3*a*b^4*e^2 + b^5*d*e))/(2*(a^7*e^7 - b^7*d^7 - 21* 
a^2*b^5*d^5*e^2 + 35*a^3*b^4*d^4*e^3 - 35*a^4*b^3*d^3*e^4 + 21*a^5*b^2*...